Resonant and Kondo tunneling through molecular magnets 
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Transport through molecular magnets is studied in the regime of strong coupling to the leads. 
We consider a resonant-tunneling model where the electron spin in a quantum dot or molecule is 
coupled to an additional local, anisotropic spin via exchange interaction. The two opposite regimes 
dominated by resonant tunneling and by Kondo transport, respectively, are considered. In the 
resonant-tunneling regime, the stationary state of the impurity spin is calculated for arbitrarily 
strong molecule-lead coupling using a master-equation approach, which treats the exchange inter- 
action perturbatively. We find that the characteristic fine structure in the differential conductance 
persists even if the hybridization energy exceeds thermal energies. Transport in the Kondo regime 
is studied within a diagrammatic approach. We show that magnetic anisotropy gives rise to the 
appearance of two Kondo peaks at nonzero bias voltages. 

PACS numbers: 73.23.Hk, 75.20.Hr, 73.63.-b, 75.50.Xx 



I. INTRODUCTION 

Over the past few years the idea of integrating the 
concepts of spintronics and molecular electronics has 
developed into a new research field dubbed molecular 
spintronicsr^ Progress has not only been stimulated by 
technological interests but has also been accompanied by 
the realization that magnetic single-molecule transistors 
exhibit various fundamental quantum phenomenal^"— 
Among many promising ideas discussed in the literature, 
particular attention has been paid to current-induced 
spin reading and writing, spin relaxation, entanglement, 
quantum computation, and Kondo correlations J^"— 

An experimental realization of spintronics devices may 
be achieved by using single-molecule magnets in com- 
bination with metallic (nonmagnetic or ferromagnetic) 
leads. For molecular-memory applications, long spin- 
relaxation times arc advantageous, which may be real- 
ized in molecules with large magnetic anisotropy, such as 
molecules based on Mni2, Fe4, and Ni4i^"— 

Controlling and detecting the molecular spin by means 
of electronic tunneling into source and drain electrodes 
poses a major challenge. While some approaches rely on 
break junctions, others are based on a scanning tunnel- 
ing microscope. In both cases, the coupling between the 
molecule and the leads can vary by several orders of mag- 
nitude, thus giving rise to strikingly different transport 
regimes .^"i^i^"^ 

In the regime of weak molecule-lead coupling, many 
experimental features such as Coulomb blockade, spin 
blockade, sequential tunneling, and cotunneling can be 
described within a master-equation or rate-equation ap- 
proach treating the electronic tunneling pcrturbativelyi^ 
However, for strong coupling, low-order perturbation the- 
ory breaks down. It is then advantageous to treat the 
electronic tunneling exactly, at the price of introducing 
approximations elsewhere. 



Recently, the Kondo effect in single-molecule magnets 
with easy-axis anisotropy has been studied by Romeike 
et alr^ Their model describes an anisotropic spin cou- 
pled to metallic electrodes by an exchange interaction, 
in the absence of a bias voltage. This differs from the 
model studied here, in which the electronic spin in the 
relevant molecular orbital is coupled to an additional lo- 
cal, anisotropic spin via an exchange interaction J, i.e., 
charge fluctuations are explicitly taken into account. In 
addition, we include a nonzero bias voltage. The presence 
of a Kondo effect for an anisotropic spin is at first glance 
surprising since it requires two approximately degenerate 
low-energy spin states connected by a term in the Hamil- 
tonian. The simplest Hamiltonian for an anisotropic spin 
S with S > I exchange-coupled to an electronic spin s 
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does not provide such a term. Using a renormalization- 
group approach, Romeike et alJ^ could show that quan- 
tum tunneling of the magnetic moment, which is de- 
scribed by higher-order anisotropy terms not included 
in Eq. ([1]), may give rise to a Kondo peak in the linear 
conductance, centered at zero bias voltage. The Kondo 
temperature is found to depend strongly on the ratio of 
the applied magnetic field and the anisotropy barrier. 
Further, Gonzalez et al^ have derived a Kondo Hamil- 
tonian of the type studied in Ref. [l^ from an electronic 
model. They have shown that a transverse magnetic field 
can induce or quench the Kondo effect. This is due to 
Berry-phase interference between different quantum tun- 
neling paths of the spin. 

Koerting et al^ consider the nonequilibrium Kondo 
effect for a double quantum dot with four leads. By re- 
moving the leads coupled to one of the dots one would ob- 
tain a model similar to ours. The main difference is that 
we include charge fluctuations on the dot which is cou- 
pled to the leads, whereas Koerting et al^ work in the 
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regime of weak tunneling and Coulomb blockade, where 
both quantum dots act as local spins. 

In the present paper, we address the question of 
spin-dependent resonant tunneling and Kondo tunneling 
through molecular magnets. As noted above, we consider 
a resonant-tunneling model where the electron spin on 
the quantum dot or molecule is coupled to an additional 
local, anisotropic spin via an exchange interaction J. We 
assume this interaction to be weak, which allows us to 
employ perturbation theory for small J . The stationary 
current through the left and right leads is identical and 
is related to the local electronic spectral function Aa{uj) 
on the molecule by the Meir-Wingreen formulai^ 

(2) 

where a is the spin, Fq is the broadening of the molecu- 
lar level due to the hybridization with lead a = L, R, to 
be defined below, and denotes the Fermi distribution 
function of lead a. The spectral function is determined 
by the imaginary part of the retarded Green's function, 
Arj{u) ~ — 2ImGJj°^(w). If transport is dominated by a 
single molecular level of energy e^, the coupling to the 
leads gives rise to a Lorcntzian form of the spectral func- 
tion, AO(c^) = T/iiuj - e^y- + T^i] with F = F^ + Tr, 
which manifests itself as a peak in the differential con- 
ductance. For single-molecule devices, the excitation of 
additional degrees of freedom due to the electronic tun- 
neling is expected to translate into additional character- 
istic features in the current. 

We consider two complementary situations. The first is 
the case of arbitrary gate and bias voltages but excluding 
the region where the Kondo contribution to the current 
is large. Within a master-equation approach treating the 
local exchange interaction perturbatively to second or- 
der, we calculate the transition rates between local-spin 
states, showing that the spin can be driven out of equi- 
librium even for strong molecule-lead hybridization. Sig- 
natures of inelastic tunneling such as the fine-structure 
splitting of the differential-conductance peaks persist in 
the regime where the hybridization energy exceeds the 
thermal energy. 

The second case concerns the regime of a large Kondo 
contribution to the differential conductance, which only 
occurs for small bias voltages on the order of \eV\ ^ 
K2{2S — 1), as we shall see. Here, transport is studied 
using a diagrammatic approach. We consider the case 
that the molecular orbital is far from resonance so that 
the resonant-tunneling contributions arc negligible. In 
addition, this allows us to obtain analytical expressions. 
We find that the magnetic anisotropy gives rise to the 
appearance of two Kondo peaks in the differential con- 
ductance at finite bias voltages ±Vc. This intrinsically 
nonequilibrium Kondo effect is quite different from the 
zero-bias peak studied by Romeike et al.^ which relies 
on higher-order anisotropics absent from our model. In 
our case, the magnetic anisotropy acts like a magnetic 



field in that it gives rise to a splitting of the Kondo peak. 
Furthermore, we find a suppression of the differential con- 
ductance with 

The paper is organized as follows. In Sec. |TT1 we intro- 
duce our model. Section IIIII considers transport within 
a master-equation approach, which allows us to study 
magnetic nonequilibrium phenomena, whereas Sec. IIVI 
considers a diagrammatic approach, which applies to the 
Kondo regime. In Sec. |V]we summarize and discuss our 
results further. Some detailed calculations are relegated 
to Appendices. 

II. MODEL 

We consider a magnetic molecule coupled to two metal- 
lic leads. Electronic tunneling through the junction is 
assumed to involve a single orbital with energy £d and 
spin s that is coupled to a local spin S via exchange in- 
teraction. The model is described by the Hamiltonian 

H ^ Hq + Hj + Hjjiiig, (3) 

where 

Hq = Ed dl^dg + Egk O'likgO-aka 

(y akcr 

+ ^ [ta al^a^a + C 4aaka) (4) 

akcr 

is the resonant-tunneling Hamiltonian, 

Hj = J s-S (5) 

with s = J2crcr' ^aif^aa' /2)(icr' IS thc cxchaugc interaction 
between the electrons in the molecular orbital and the 
local spin S, and 

iJmag = -K2{S'y (6) 

describes the easy-axis magnetic anisotropy of the local 
spin. We choose the z axis as the easy axis. Here, 
creates an electron with spin a and energy Sd on the 
molecule, while a^j^.^ creates an electron with energy e^k, 
wave vector k, and spin a in lead a, which is considered a 
noninteracting electron gas. The vector cr = {a^, (Jy,(Jz) 
denotes the Pauli matrices. In break junctions produced 
by electromigration, the onsite energy Ed can be tuned 
by applying a gate voltagci^"— 

III. MASTER EQUATION FOR THE SPIN 

The presence of strong coupling between the molecule 
and the leads prevents us from treating the hybridiza- 
tion term in Eq. (|4]) perturbatively. However, since the 
Hamiltonian becomes bilinear in the limit of vanishing 
exchange coupling, J = 0, our strategy is to diagonalize 
Hq + H-caag exactly while treating Hj as a perturbation 
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up to second order. This approach ahows us to study 
the nonequihbrium dynamics of the molecular spin at 
finite bias voltages for arbitrary molecule-lead coupling 
strengths, provided that Kondo correlations do not lead 
to a diverging contribution from higher-order terms in 
the expansion. 

We start by rewriting Hq in terms of new opera- 



tors 



35.36 



where 



a'k' 

d^a — ^ ^ ^^QkCak(T; 
ak 



and 



Eak - Ca'k' 



Eak - £d- J2a' 



(7) 

(8) 
(9) 

(10) 
(11) 



For simplicity we assume real tunneling amplitudes ta- 
in terms of the new operators, the exchange interaction 
assumes the form 

Hj ^ J '^ak'^a'k' clk^^y-Ca'k'<T' ' S. (12) 

aa'kk'crcr' 

The time evolution of the density matrix p of the full 
system is described by the von Neumann equation, p = 
— {i/h)[H, p]. The degrees of freedom of the local spin 
are described by the reduced density matrix 



pj = Trci p, 



(13) 



which is obtained by tracing out all electronic degrees of 
freedom. Assuming that the large electronic subsystem, 
which acts as a spin reservoir, is weakly perturbed by 
the exchange coupling, we replace the full density matrix 
by the direct product p ~ p,/ pc\. We need a further 
approximation for the electronic density matrix pei- We 
assume that different chemical potentials p^ (p-fi) are 
imposed for the left (right) lead far from the junction. 
It would thus be natural to assume Fermi distributions 
fa{uj) = f{uj ~ Pa) for the physical electrons created by 
^Iko-- However, we need to make a reasonable assump- 
tion on the transformed c fermions appearing in Eqs. ([7]) 
and (|T2|) . Since c^^j.^. (cj^j.^.) creates an electron in a state 
with vanishing probability density far from the junction 
in the right (left) lead, we assume the occupation num- 
bers of these states to be described by /q(w). 

Making use of the Markov approximation that pj 
changes slowly on the time scale of electronic relaxation, 
we obtain 



1 /■* 



Pel]' 



Here, operators 0(t) with an explicit time argument, 
including pj{t), are in the interaction picture, 0(t) = 
e^{Ha+Hr^^^)t/h Q ^-^(Ho+H-^.^^)t/h _ ^^^^ ^YiaX sccoud-or- 

der perturbation theory in the exchange coupling gives 
the first non-vanishing correction to the conductance, 
since the expectation value (S) and thus all first-order 
terms vanish exactly due to symmetry. 

We are interested in the stationary state. The station- 
ary density matrix pj has to be diagonal in the basis of 
eigenstates \m) of S*^, since the full Hamiltonian H is 
invariant under rotation about the 2-axis in spin space. 
Inserting Eq. into Eq. ([H)) we thus obtain a Pauli 
master equation, also called rate equations, of the form 

— Pm {Rm^m+1 + Rm^m-l) ~ (15) 

for the occupation probabilities P,„ of spin states |m) in 
the stationary state. The transition rates read 



Rn 



2±1 = i 115*^1771) 



|2 

2TTh 



a.a.' 



dw |i>a(w)|^ |z>a.(w - [±2m + l]K2}f 

X [1 - faiio)] fa'iuj - [±2m+ 1]A'2). (16) 
The spectral functions are given by 



|i>Q(w) 



(L.-ed)2+rV4 



(17) 



with r = Fi -I- Ffl and F„ = 27ri^£)^. The densities of 
states for the leads. Da, are taken as constants. Com- 
pared to Eq. (|lip we have approximated the self-energy 
part of i^ak by a constant and absorbed a factor 2-KDa- 

At zero temperature, the integrals can be evaluated 
analytically. In the limit of large bias voltages, the rates 
approach the constant value 



i?. 



Ll R 



■rn— »?7i±l 



X |(to±1|5±|?7i)|2. (18) 

On the other hand, at zero bias only the rates involving 
the absorption of energy arc finite, whereas the emission 
rates vanish. 

Solving Eq. (fTS)) allows us to compute the differential 
conductance of the molecular junction. The current op- 
erator of lead a reads 

e 

k 

.e 

'ft. 



kcr 

= 4 II (^"<'k'^a"k"C^kvCc."k"<. - h.. 



kcr Q'k'Q:"k" 



(19) 



(14) 



In order to compute the spin-dependent contribution to 
the expectation value (/„} = Tr/^p = Tr Ia{t)p{t) of 
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the total current, we use the iterative solution of the von 
Neumann equation, 

t t' 

P(t) = -4 / I ^Hj{t'). [Hj{t"), p{t")]] . 

A term containing p(—oo) has dropped out here since it 
is linear in Hj and thus vanishes upon taking the trace. 
Making use of the Markov approximation we find 



dt' / dt" 

: J — oo 

^Tv[[Ut),Hj{t')],H.j{t")]p{t) (21) 



for the second-order term. Carrying out the time inte- 
grals and evaluating the spin sums as explained in Ap- 
pendix |Al we obtain 



{II) 



(2) 



2TTh 4 



1 



|(m - l\S-\m)flc,a'a"{[-2ni + IjA'a) 



+ 2 K"^ 



+ \{m\S'\m)fl, 



|m)|2/„ 



,{[2m + l]K2) 



'(0) 



(22) 



with 



Iaa'a"{E)= J duj 

X {|;>„(l^)P|£>,,(c^-£;)|2 [1 - f^{uj)] {lu - E) 

- |z>a(w)P |z>a'(w + E)\^ /„(w) [1 - /,.(l^ + 

(23) 

Equations and ([25)1 give the first non-vanishing cor- 
rection to the zero-order current which is ob- 
tained from the Meir-Wingreen formula [Eq. by in- 
serting the spectral function of the unperturbed system, 
A°{uj) = r/[(e<j - w)2 + rV4]. Note that Eq. @ with 
= is recovered by inserting the equilibrium den- 
sity matrix p° and the current operator from Eq. 
into = Tr/ipO. 

A simple interpretation of Eq. ([^^ is possible for the 
special case of a local spin of length 5 = 1/2, for which 
the magnetic anisotropy K2 is irrelevant and can be set 
to zero. For this case we obtain 



(2) ^ J'^SjS - 



duj 



1) TrT 



1 3 



(c^-ed)2 + rV4 



[/i(c.)-/fl(c.)]. (24) 



Here, the third power of the spectral function ap- 
pears, since the current operator and the two exchange- 
interaction operators in Eq. (|2ip are each bilinear in 
fcrmionic operators. 
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FIG. 1: (Color online) (a) Current-voltage characteristics for 
different hybridization energies, F = K2/2Q, F = 7^2/10, and 
F = K2- The inset shows a closeup of the fine structure 
at positive bias, (b) Magnetic transition rates 7?2-»i, Ri^o 
and occupation probabilities P,„ as functions of bias V for 
F = A'2/20. We assume symmetric couplings to the leads, 
Tl = F_R, and symmetric capacitances, pL ~ eV/2, pR — 
—eV/2, a local molecular spin of length S = 2, Sd ~ 4_fL2, 
and zero temperature. Further, we set J = F/5. Currents are 
given in units of (2e/fi,)Fi:,Fjj/F. Rates are given in units of 
their maximum values, cf. Eq. (I18|l . 



If the magnetic anisotropy is large compared to the 
hybridization energy, A'2 ^ F, the general expression for 
the current in Eq. (|22p simplifies to 



^^>*'^= E(ri-^^"[ri + ri?])W'(o) 



^P,„|(TO|5^|m}| 



(25) 



since the integrals /q,q/(j// ([±2m + 1]K2) are negligible 
compared to /qq'q"(0). Assuming symmetric capaci- 
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FIG. 2: (Color online) (a) Current-voltage characteristics for 
different hybridization energies, F = 7^2/20, F = K2/IO, and 
F — K2- The inset shows a closeup of the fine structure 
at positive bias, (b) Magnetic transition rates i?2->i, Ri^o 
and occupation probabilities P,n as functions of bias V for 
F = _fL2/20. We assume strongly asymmetric couplings to 
the leads, Fl <^ F_b, and strongly asymmetric capacitances, 
fj,L = eV, fiR = 0, a local molecular spin of length 5 = 2, 
Ed = 4X2, and zero temperature. Further, we set J = F/5. 
Currents are given in units of (2e/ft)FiF_R/F. Rates are given 
in units of their maximum values, cf. Eq. psp . 



tances, one finds, in tlie limit of large bias voltages, 



2e TlTr 

h r ' 

2e TlTr ^S{S 
h T 4 



p2 ' 



(26) 
(27) 



for the zero-order and second-order contribution, respec- 
tively. Note that the ferromagnetic or antiferromagnetie 
sign of J does not affect the results in the present ap- 
proximation. 

We first consider the situation of symmetric molecule- 
lead couplings and capacitances, i.e., \tL\ = Fl = 
Tr, and /ii = eV/2, fin = —eV/2. Figure [TJa) shows 
the current-voltage characteristics up to second order in 
J for the case of a local spin of length S = 2. The 



characteristic fine structure of the current step at the 
Coulomb-blockade threshold is due to the second-order 
contribution, whereas the main step is mostly 

coming from The fine structure persists as long 

as the hybridization energy F remains small compared to 
the magnetic anisotropy K2- Note that the broadening 
of the steps is due to F > 0, and not to the temper- 
ature, for which we assume T T. For bias voltages 
below \eV\ — 2ed, the current and all magnetic excita- 
tions are thermally suppressed. However, as soon as the 
chemical potential of one lead aligns with the resonance 
of the molecule, the current increases to its maximum 
value. The current-induced magnetic transitions become 
energetically accessible at the same time, as shown in 
Fig.[TJb), resulting in nonequilibrium probabilities Pm of 
the different spin states. In the limit of large bias voltages 
all spin states are equally occupied, P,„ = 1/(25* + 1). 

Interestingly, the presence of magnetic anisotropy leads 
to negative differential conductance in the vicinity of 
\eV\ ~ 2ed- The underlying mechanism shall be ex- 
plained briefly. According to Eq. ([25|) . the magnetic 
states with maximum quantum numbers m = ±5* domi- 
nate the current since the current is proportional to the 
average '^TnPm\{'m\S'^\m)\^ . Each decrease in P±s thus 
causes a decrease in the current. Therefore, the spin- 
dependent contribution to the current is large at low 
bias voltages, where P±2 = 1/2 and X^m ^"i"^^ " 4, 
and small at high bias voltages, where P±2 ~ 1/5 and 
^^P„m2 = 2. 

We next consider the situation of strongly asymmet- 
ric molecule-lead couplings and capacitances, i.e., I^ij ^ 
ifll, Tl Tr, and ^ml — eV, fin ~ 0. Note that 
this regime is naturally realized in many experimental 
setups, whereas perfectly symmetric couplings are in gen- 
eral more difficult to achieve. The current- voltage curves 
and the bias dependence of the magnetic excitation rates 
are shown in Fig. [51 Due to the asymmetric coupling 
and Ed > 0, the current is suppressed for negative bias 
voltages. However, the characteristic steps correspond- 
ing to excitations of the molecular spin reappear at pos- 
itive bias. Only their abscissas are reduced by a factor 
of 2, since the chemical potential of the left lead is now 
fiL = eV instead of fiL = eV/2. The device thus acts as 
a rectifier. Note the small ohmic contribution with con- 
stant slope for large coupling F in Fig. [DJa). We return 
to this point below. 

We finally turn to the full bias and gate-voltage de- 
pendence of the current. In both the symmetric and 
the asymmetric case, selection rules for the spin require 
changes in the magnetic quantum number by Am = 
or Am = ±1, where Am = corresponds to elastic and 
Am = ±1 to inelastic scattering events, cf. Fig. [31[a). In- 
elastic scattering processes appear as additional steps in 
the current and give rise to the magnetic fine structure 
in the two dimensional density plots of the second-order 
contribution to the current as a function of bias and gate 
voltages shown in Figs. [2Jb), (c). We can now under- 
stand the origin of the weak ohmic conduction seen in 
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FIG. 3: (Color online) (a) Level scheme showing all spin tran- 
sitions to order J^. (b),(c) Two-dimensional density plots of 
the absolute value of the second-order current contribution 
(7l)'^' as a function of bias voltage V and gate potential 
for r = Kij'Xi- Ed is controlled by the gate voltage. We 
choose the same parameters as in Fig. [T] Bright (dark) colors 
correspond to high (low) currents. In (b) we assume sym- 
metric couplings, hl ~ eV/2, fiR = —eV/2, Fl = Tr, while 
in (c) we assume asymmetric couplings, fiL ~ eV, fiR = 0, 
Tr > Pi. 



sary to resum the divergences to all orders in J42rii 

The total current through the molecule is related to the 
local electronic spectral function Aa{uj) = — 2ImG'J'j^*(w) 
by the Meir-Wingreen formula, Eq. © , where G™^/ (w) = 
/ d(t - t') e"^(*-*')G"JJ,', {t - t') denotes the Fourier trans- 
form of the retarded single-particle Green's function 



(28) 



Making use of the transformation defined in Eqs. ([5|)- pT|) 
requires to compute the finite-temperature time-ordered 
Green's function 

Gaa'kk'aa'iT,T') = - (Tr CaWa {t)cI,-i^, ^, {t )) . (29) 

Our strategy is to expand Gaa'kk'aa' in powers of J. 

In order to obtain the current from the Meir-Wingreen 
formula, we need the imaginary part of the electronic 
Green's function, 

Im^G-*(..)=Im Yl '^ckK'U'Gal'kU'.A^), (30) 

where G^aa'kk'aai'^) denotes the retarded Green's func- 
tion. All non-vanishing diagrams of the local Green's 
function up to third order in J are shown in Fig. |4j fol- 
lowing the notation of Ref. [H. We again consider the 
case of strongly asymmetric couplings, i.e., [t^l ^ \tji\, 
^ ^R, and — eV, fiR ~ 0. As we shall see, 
Eq. pop is then dominated by the contribution from the 
right electrode, a = a' — R. This allows us to describe 
the molecular degrees of freedom by a thermal equilib- 
rium distribution function that is independent of the ap- 
plied bias voltage and to obtain an analytical expression. 



Fig. [2ja) for F — K2. What we are seeing is the tail of 
the current step at = in Fig. ^c), which is consid- 
erably broadened for F = K2- Note that we here have 
ed = 4:K2 = 4r, i.e., we are only 4F away from the step. 
Since this distance does not depend on the bias voltage, 
the conductivity is essentially constant, leading to ohmic 
behavior. 

IV. KONDO TRANSPORT 

Second-order perturbation theory in the exchange in- 
teraction J fails, even for small J, if the prefactors of 
higher-order terms diverge. This is the case in the Kondo 
regime. Logarithmic divergences of the conductance first 
appear in terms of third order in JiiSr— (In this section, 
we assume antiferromagnetic exchange, J > 0.) Studying 
the emergence of Kondo correlations thus requires to go 
beyond the second-order master-equation approach dis- 
cussed in Sec, mil For sufficiently small J and sufficiently 
large thermal energies, the conductance is dominated by 
the third-order contribution, which we calculate in this 
section. At lower temperatures, it would become neces- 
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FIG. 4: Non- vanishing diagrams of the impurity Green's func- 
tion up to third order in J, following Ref. [3^ . Diagrams in- 
cluding fermion loops vanish exactly and are not shown. Spin 
averages are denoted by dotted lines. 

The evaluation of Eq. ([50)1 is shown in Appendix |B] 
We obtain 

Im V G:,<^(w) ~ — , , + ( -^^-^^—] ImI]'°Vcj) 

(31) 

with 
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T vrct/ N ^^72 / nV^P 1 - f {oJ + E„i - El) j ^ 

ImS (w) = - -J iyo{ed)2_^P,n ^ _ ^ )^ni2_^ 

mnl -'^ ' { i 



|(m|5»|^ 



ijk 



In 



+ ln 



(32) 



where 



0.6rG(e2//;) 



t'o(ed) = 



r/27r 



r2/4' 



(33) 



In the derivation we have assumed the molecular level 
to be far from resonance, i.e., \ed\ to be large compared 
to K2S, T, and F, but still small compared to the band 
width X of the leads. Details are discussed in Appendix 
IbI We have also assumed |w| ^ le^l, the significance 
of which will become clear in the following step. Under 
these conditions, the resonant-tunneling contribution to 
the differential conductance, which we have studied in 
Sec. mil is negligible compared to the Kondo contribu- 
tion. 

In the low-temperaturc limit T <^T, derivatives of the 
Fermi functions with respect to the bias voltage become 
delta functions and the differential conductance simplifies 
to 



G 



2Tih F 



rV4 



F2/4 



(34) 



Note that the argument lo of 17'^^ [lo) is eV . The as- 
sumption |a;| <^ \ed\ made above thus corresponds to \ed\ 
also being large compared to the bias, \eV\. The spec- 
tral function has logarithmic divergences for T — at 
the transition energies of the molecule, Em — Em cor- 
responding to virtual transitions between two magnetic 
states m) and |n). One recovers the prefactor SttJ'^Do/S, 
see Ref. |3^, for the case of spin 5 = 1/2 and the (then 
irrelevant) anisotropy set to K2 = 0. 

Numerical results for nonzero temperatures are shown 
in Fig. [5] The differential conductance G diverges loga- 
rithmically for r ^ at critical bias voltages V — zkVc 
with eVc = Es-i —Es = K2{2S —1) since the emergence 
of Kondo correlations requires the bias voltage to ex- 
ceed the energy of the transition from the ground states, 
TO = zLS, to the first excited states, to = zt{S — 1). 
Note that G{V) is symmetric for positive and negative 
bias, in spite of the highly asymmetric coupling since it 
is probing the electronic spectral function. The situation 
is quite different from the case considered by Romeike 
et al.^ which concerns a zero-bias peak resulting from 
quantum tunneling between the two states \S) and |— S"). 



r=o.i 




-10 
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FIG. 5: (Color online) Differential conductance G in units of 
e^//i for different thermal energies T (in units of K2) obtained 
from Eq. (|34|) as a function of bias voltage V in units of K2 /e. 
Here we assume a local molecular spin of length 5 = 2 and 
choose Jvo{ed) = I, x = 100i4:2, and Tr = 100 Fl. Note 
that the parameters F, Ed, and J leave the curves for G (in 
arbitrary units) invariant except for changing the constant 
offset. 



In our case, the splitting of the Kondo peak as a con- 
sequence of magnetic anisotropy is more similar to the 
situation of a quantum dot in an external magnetic field 
with Zeeman energy B, where a Zeeman splitting of the 
energy levels leads to the occurrence of two conductance 
peaks at eV ~ ±B in the Kondo regime4^ At higher 
temperatures, T ^ K2, the two Kondo peaks merge into 
a single peak centered at zero bias due to the thermal 
excitation of spin states with higher energy. 

We now turn to the Kondo temperature Tk- Poor 
man's scaling for the equilibrium case results in Tk = 
because the matrix elements of between the two 
degenerate ground states of the local spin vanish for our 
model if S" > l/2iiS Since a Kondo effect evidently does 
occur at nonzero bias, this result is clearly not sufficient. 
A rough estimate of the Kondo temperature Tk can be 
obtained as the temperature for which the second-order 
and third-order terms become equal in Eq. ([5^ . We find 
Tk ^ eyi-p[—\/avQ{ed)J], where a is a number of the 
order of unity. In the limit K2 0, where the two peaks 
in Fig. [S] would merge, we recover the result a = 2 for an 
isotropic spin. 
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Since we focus on the case of strongly asymmetric cou- 
plings, where the molecular degrees of freedom are in 
equilibrium with one of the two leads, the logarithmic 
divergences are cut off by temperature or the applied 
bias voltage, respectively, in our perturbative approach, 
see Eq. ([32]) . The divergence for T — !■ is unphysi- 
cal and would likely be removed by a resummation of 
higher-order terms. By analogy to Ref.l45l we conjecture 
that the divergence is ultimately cut off by a voltage- 
dependent spin-relaxation rate. 

While we have so far discussed the dependence on the 
bias voltage, see Fig. [51 we now turn to the gate voltage. 
The gate voltage shifts the on-site energy Sd and thus 
enters the expression for the current through the square 
of the spectral function r/[£^ -I- r^/4] and the square 
of the factor ed/ie^ + r^/4]. In particular, we obtain a 
suppression of G cx in the limit of strong detuning. 



First, it would be interesting to include a local Coulomb 
interaction U between the electrons on the molecule. 
However, due to the large hybridization there are no 
states with large probability on the dot and the effect of 
U is expected to be relatively weak. We expect that for 
very large U an equilibrium Kondo resonance could oc- 
cur as a zero-bias peak in the differential conductance in 
addition to the nonequilibrium Kondo effect described in 
this paper. Second, the presence of an external magnetic 
field might lead to an interesting interplay with the split- 
ting of the Kondo peaks due to the magnetic anisotropy. 
Finally, it would be desirable to combine the two cases 
studied here and to analyze the Kondo effect in mag- 
netic molecules in the resonant-tunneling regime, where 
resonant-tunneling contributions to the conductance are 
not negligible and the spin is driven out of equilibrium 
by the current. 



V. CONCLUSIONS 

We have studied the spin-dependent electronic trans- 
port through magnetic molecules for strong coupling to 
the leads. Our discussion has focused on two complemen- 
tary regimes. 

For the first regime, we have presented a description 
of transport in terms of a master equation that keeps the 
electronic tunneling exactly, holds for arbitrary bias and 
gate voltages, and treats the local exchange interaction J 
perturbatively at second order. This approach is thus ap- 
plicable for small J. We have derived the bias-dependent 
magnetic transition rates showing that the tunneling cur- 
rent can be used to drive the molecular spin out of equi- 
librium. Further, we have shown that the characteristic 
fine structure of the differential-conductance peaks per- 
sists for strong molecule-lead coupling, where the broad- 
ening of the peaks is determined by the hybridization 
energies. 

The perturbative expansion in J fails if Kondo corre- 
lations contribute significantly to the transport. In this 
case, prefactors of the third- and higher-order terms in 
J diverge for T — > 0. The Kondo correlations can be- 
come important for small bias voltages on the order of 
\eV\ ^ K2{2S — I). Here, transport is described by the 
Meir-Wingreen formula in combination with a diagram- 
matic calculation of the local electronic spectral function 
of the molecule. We have assumed the molecular level to 
be far from resonance, which on the one hand makes sure 
that the resonant-tunneling contributions to the conduc- 
tance are small and which on the other allows us to obtain 
analytical results. We have shown that Kondo peaks ap- 
pear at finite bias voltages proportional to the anisotropy 
energy of the molecular spin. 

Our results leave several avenues for future research. 



Acknowledgments 

We would like to thank A. Donabidowicz-Kolkowska, 
D. R. Reichman, and A. J. Millis for useful discussions. 
Financial support by the Deutsche Forschungsgemein- 
schaft is gratefully acknowledged. 



Appendix A: Calculation of the current 

In this appendix we give details on the derivation of 
Eqs. (HH) and i^. We start from Eq. 



" J —oo J —I. 



yiTv[[Ut),Hj{t')],H.j{t")]p{t). (Al) 

Inserting the expressions for the current operator /„, 
Eq. (fT9|) . and for the exchange interaction i/j, Eq. ()12|) . 
we find 



' h Ah? 



j dt' j dfTiY, -^--i^-i 



kcr 123456 



cl(t)c2(t),4(i')c4(t')cr.3^4 ■S(t') 



4(f")c6(i")Ta..e ' S(i") 



(A2) 



where we have assumed to be real. Here, the short- 
hand notation j ~ 1, 2, 3, 4, 5, 6 stands for (aj, kj, aj). 

Introducing t = t — t' and t' = t' — t" and assuming a 
product state gives 
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, r2 poo poo 

•^0 -^0 k 123456 

X { [623545SM4 (1 - /e) e'(-'-'-^)-/^e'(-<>--^y/^ - SuS25S3e.f2 (1 - .A) (1 - /e) e'(-^-'')-/^e'(-'-''y 

- 5u52M2fA (1 - h) e^i-^-'-^-l^e^'^-'-^^'l'' + ,523<545^i6 (1 ^ h) /4 (1 - k) e'(^=-^^)-/''e^('«-'^)-'/?>] 
X Trj2S(0) •S(-r')pj 

- [<523<^45<5l6/2 (1 - /4) /ee'(^e-^=)-/''e'(^«-^*)-'/'i _ S^^^^^b^^ (l _ J^) (l _ J^) /ge'(e4-e.)r/;.g,(.,-.,)r'/?. 

- <5i4<^25<^36 (1 - h) f^fee'^---^)^/^e'^-<>---y/^ + 623645616 (1 - /2) (1 - /4) /6e'(^^-'^)-/''e»(-«--^)-'/?>] 

X Trj2S(0) •S(r')pj}. (A3) 
This result can be rewritten as 

•'^ ''^ k 123456 

X { [523^45-5x6/4 (1 - h) e'(^«-^^)-/ V^^^"^^)^'/" - <Si4<525<536/2 (1 - k) e*(^4-^=)-/ V^^^"^^)^'/^] Tr j 2 S(0) ■ S(-t') pj 

- [<523<545<5i6 (1 - h) /6e*(^«-^^)^/''e'(^«-^*)^'/'' - (5i4<525<536 (1 - /2) /6e'^(^^-^=)^/''e'(^«-'^)^'/''] Trj 2 S(0) ■ S(r') pj}. 

(A4) 

I 

The sums over spin indices are simphfied by making use Noting that Eqs. (jlOp and pip imply 
of the identities 

a „ c ocjc ^ 6La^+mDLtLVc,^{ei) (A7) 

2_^a-cra'-^l CTcrV ■ »2 = ^^1 ' O2, ^ 

aa' 

^ctctct'-Si ctct', _cr • S2 = 0. (A5) and 

era' 

* r 2 

In the coefhcients v^i^ in Eq. ((TT]), we approximate the Va{e) — Va{e) ~ ^~ ' 

self-energy part by a constant, as we did in Sec. IIIIl 

^ we arrive at the following expression for the tunneling 

i^Qk = t^a(eak) ^ — TTTfK- current: 

Cak - Ed - «i /2 

I 



{^l)^"^^ = '^T—TTT I / E [i^Lai - i'^DLtLvl)v2vlViVlvQ - vl{6La., + iT^DLtLV2)vlv4vlvQ\ 

^ -^0 -^0 123456 

X { [523SA55l6h (1 - k) e'(^«-^=)-/ V(^«-^^)^'/'"' - <5l4<525536/2 (1 " /6) e^('4-..)r/ V(^6-.2)r7;.] Tr, 2 S(0) • S(-r') 

- [<523<545'5i6 (1 - h) /6e*(^«-^^)^/''e*(^«-^^)^'/?» _ 514^25536 (1 _ f.,) f^e'i-^---)r/^e'i-e-e.)r' /n^ 2 S(0) • S(r') pj}. 

(A9) 

I 

Since we have assumed pj to be diagonal in the stationary Appendix B: Calculation of the impurity Green's 
state, we finally obtain Eqs. ([^ and ([221) • function 



In order to use the Meir-Wingreen formula for the con- 
ductance, we have to compute the imaginary part of the 
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Green's function in Eq. ((30)) . We consider the situation 
of strongly asymmetric molecule-lead couplings and ca- 
pacitances, i.e., It^l Fi ^ Tji, and /il — eF, 
A^i?. - 0. 

Since Wick's theorem does not apply to spin operators, 
averages of products of spin operators do not factorize 



into averages of pairs. We follow Ref. |38| in evaluating 
the spin averages. Expanding the electronic Matsubara- 
Grcen's function in powers of J and organizing the ex- 
pansion in terms of topologically distinct diagrams, one 
obtains^^ 



a„o'kkV.'(r,r') = ^ — r dr, ■ ■ ■ dr^ E (t. [^^^ (n) • • • 5'" (r„)] ) 



X (Tr 



(Bl) 



where /3 = 1/T denotes the inverse thermal energy. For 
convenience, we have defined i?cr = X^ak ^akCaka- All 
non- vanishing diagrams up to third order in J are shown 
in Fig. m The linear term vanishes, since (S) = 0. 
Diagrams with fermion loops are zero for the following 
reasons)^ a loop with a single fermion line results in 
taking the trace of the Pauli matrix in the vertex, which 
yields zero. A loop with two fermion lines appearing in 
the third-order diagrams gives rise to a trace over two 
Pauli matrices, TrcrV-' = 2(5,,,. The resulting spin aver- 
age {Tr[S'^{Ti)S'^{T2)S'^{T3)])o, with at least two of ij, 
i2, and 13 equal, vanishes. 

Splitting off the zero-order term, the Green's function 
in Eq. ()Bip can be written as^ 



Gaa'kk'acr'{T,T') = Gakai'''^ '''' )^aa' Skk' S^a' 

+ dn dT2gli^^{T,Ti)i:aa'kk'aa'{Tl,T2) 

Jo Jo 

xekV'(^2,r'), (B2) 

where the unperturbed Matsubara-Grecn's function in 
the imaginary-time domain is given by 

Gli^Ar, r') = - [e{T r') /K^)] 

(B3) 

with cJak = Ectk — fJ-a- In the frequency domain we have 



IUJ„ — t^Qk 



(B4) 



where iuj„ is a fermionic Matsubara frequency. Note that 
we are only interested in the spin trace of the self-energy, 
So- ^aa'kk'o-cr' ^hich cntcrs in the Mcir-Wingrcen for- 
mula. 

The second-order term of the self-energy yields 

cr ctikiCTi m,n 

X J^*ki^a'k'|i^aikiP 

X ^ |(m|5»pe(^'"-^")(— ^)/^P„„ 



(B5) 



where we have used that TrcrV-' = 2Sij. Here, = 
—K2'm^ denotes the magnetic anisotropy energy in the 
spin state |m) with occupation probability Pm- The spin 
averages in Eq. (jBl[) arc to be evaluated for the unper- 
turbed Hamiltonian Hynagr^ leading to P,„ oc e~^^"^. 
We restrict ourselves to the off-resonance situation, i.e., 
the dark region in Fig. ^c) with negligible resonant- 
tunneling differential conductance, where the spin essen- 
tially remains in equilibrium. This is certainly satisfied 
if \ed\ is large compared to the energy scales relevant for 
the Kondo contributions, K2S and T. 

Equation (|B5p contains a sum over leads, ai ~ L, R, 
and a factor of t^^ under the sum. Since we have assumed 
strongly asymmetric couplings, <^ |tij|, the sum is 
dominated by the contribution from the right lead, ai = 
R. Dropping the term with ai = L, we note that the 
Green's function Q^x^rri'^^T') in Eq. I|B3|) only contains 
the Fermi distribution function for the right lead, which 
is /H(e/?k) = /(cflk) = l/(e'^<=«'' + 1), since ^ir = 0. 
Importantly, the resulting expression is independent of 
the bias voltage. 

Furthermore, we see that Eq. (|B5[) contains a factor 
tata'- From Eq. ([50)) we obtain the same factor so that 
the contribution from leads a, a' is proportional to t^i^/. 
Since we have assumed |ti| ^ |tfl|, we can neglect all 
contributions except for a ~ a' = R. We will keep only 
these contributions from now on. 

Taking the Fourier transform of Eq. (|B2[) and perform- 
ing the analytic continuation, we obtain the retarded 
Green's function 



1 



p- 



1 



— in 5{uj — e/fk 



ei?,k' 



(B6) 



where P denotes the principal value. We assume le^l to 
be large not only compared to K2S and T but also to 
F. One can then show that the delta-function terms are 
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negligible compared to the principal value terms. Includ- 
ing the factors of t^lfj,., vru., we obtain expressions of the 
form 



k 



1 



R 



uJ-eRk r (erf-w)2 + r2/4 



(£d-^)2 + r2/4- 



(B7) 



For the imaginary part of the Green's function in Eq. (j30p 

we then only require the imaginary part of ^flflkkVo- 
in Eq. (|B6|) . Taking the imaginary part of the Fourier 
transform of Eq. (jB5[) we obtain 



ttJ^Do 



X'^flk^-Rk' \vr{u; + Em - E„ 



,2 1~/(W + S„ 



E„ 



1 - /(c^) 



(B8) 



where we assume constant densities of states for the 
leads, Dq = = Dr, and an energy band ranging 



from —X to X, where x is the largest energy scale in our 
model. 

The third-order term gives 



II^fl]?kkw(^l'^2) ~ ^3 



J3 



<y ki£Ti,k20-2 



/•/3 

ki?kiPkflk2p / '^7-3 0Hki^i('^l''^3)^V,<T2(^3,T2) 

Jo 



ijfc 



[5'(Ti)5^"(T3)5'^(r2)] 



CJ"i(T2 <72Cr 



2 



(B9) 



Here, the average involving spin operators depends on 
the time arguments Ti, and r3, since i, j and fc can 
be different. However, since the self-energy only de- 
pends on the differences t\ — T2 and T3 — ri, we may 
set T2 = and distinguish the two possibilities ti > T3 
and T3 > Ti. Using that Tr[<TV-'<T'''] = 2ieijk, inserting 
^V(n,r3) = -Kri-r3)-/(eflk)]e-^«''(^i-^-^) and 
^V(^3,0) = - [0(t3) - /(e/jk)] e-^^-^^S and evaluating 
the integral over , we obtain for < ri < /3 



J 



E4Vw(^i.o) = E '^flk^wkflk,ni^iik2i'E^»^'=E("^i^»("i^-''io(^i5'i"^)^" 

O" kik2 mn/ 



ei?ki — eflk2 + En ~ 




1 - /(Cflks) 




EiJki — Eflka + - 




/(ej?ki) 




ei?ki — eHk2 + E,n ' 


- En 








■[l-/(ei?k2)]e-^-''^"^/V^--^')"^/'' 
^ [1 - /(ei^kj] e-^-i-i/''e(^"-^''^^/4. (BIO) 

Eflki — e^ka + J^rn — J^n ) 



With {m\S'\n)*{n\S^\l)*{l\S''\m)* = -{m\S'\7j){n\S^\l){l\S''\m) under the sum over i, j, k, Eq. (|BT0)) simplifies to 
E^flflkkVa('^i'O) = E '^ki'i?kni^flkilVi?k2rE'^y'=E^'^l'^'l")*^"l'^^l')^^l'^''l'^)-^™ 

kik2 'ijA; mnl 



(Bll) 



Computing the Fourier transform yields 

E^fi«kk'..(^) = E^''^-fcE(™i^»("i^'io(^i^'i"^)^™ E ^k^i?k'i^flk.pkHk2p 

<7 'i ?fc mn^ ki ko 



[ l-/(ei?kj 1 l-./(£flk2) 

I ei?ki — eflk2 + En — El uj - e^jka + Em — El +iS 1 - f{€Rj^.^ - Em + -Bi) 

/(cflki) 1 1 - ./(£i?k2) 

ei?.ki - eflk2 + Em - En us - eru.^ + Em - El +iS 1 - /(efik2 - Em + Ei] 



(B12) 
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The sum over k2 can be evaluated to give 



a 

x\vr{uj + E„, - Ei)\^ 



J2^^jkY.HS'\n){n\Sm){l\S''\m)P^J2\ 



ijk rani 

l-/(ei?kj l~ f{w + E,^~Ei) 



ki 



I- f{LO + E,r,-Ei) 



Cflki — W + i?„ — i?„ 



1 - /(^) 



e/v^ki — uj + El — En 



1 - /(^) 



r 



(B13) 



Finally, the sum over ki leads to Eq. ([32|) for the self- n and only keep the terms that diverge at cj = En — E„ 
energy. Here we assume x 3> |ed| ^ ijJ,En for all states and low temperatures. 
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